Spin canting in a Dy-based Single-Chain Magnet 
with dominant next-nearest neighbor antiferromagnetic interactions 
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We investigate theoretically and experimentally the static magnetic properties of single crystals of 
the molecular-based Single-Chain Magnet (SCM) of formula [Dy(hfac)3NIT(C6H40Ph)] oc compris- 
ing alternating Dy''''"and organic radicals. The magnetic molar susceptibility xm displays a strong 
angular variation for sample rotations around two directions perpendicular to the chain axis. A 
peculiar inversion between maxima and minima in the angular dependence of xm occurs on increas- 
ing temperature. Using information regarding the monomeric building block as well as an ab initio 
estimation of the magnetic anisotropy of the Dy^+ ion, this "anisotropy inversion" phenomenon can 
be assigned to weak one-dimensional ferromagnetism along the chain axis. This indicates that an- 
tiferromagnetic next-nearest-neighbor interactions between Dy^+ ions dominate, despite the large 
Dy-Dy separation, over the nearest-neighbor interactions between the radicals and the Dy^"*" ions. 
Measurements of the field dependence of the magnetization, both along and perpendicularly to 
the chain, and of the angular dependence of Xm in a strong magnetic field confirm such an inter- 
pretation. Transfer matrix simulations of the experimental measurements are performed using a 
classical one-dimensional spin model with antiferromagnetic Heisenberg exchange interaction and 
non-coUinear uniaxial single- ion anisotropics favoring a canted antiferromagnetic spin arrangement, 
with a net magnetic moment along the chain axis. The fine agreement obtained with experimental 
data provides estimates of the Hamiltonian parameters, essential for further study of the dynamics 
of rare-earths based molecular chains. 

PACS numbers: 75.50.Xx,75.30.Gw,75.30.Cr,75.10.Pq 



I. INTRODUCTION 



For many years, one-dimensional (ID) magnetic sys- 
tems have been intensively studied owing to their sim- 
plicity. A number of sophisticated theoretical predic- 
tions have been experimentally verified as, for exam- 
ple, Haldane's conjecture for Heisenberg antiferromag- 
nets with integer spin^, Villain's mode for quasi-Ising 
s = 1/2 antiferromagnets^, the contribution of solitons 
to thermodynamic properties for easy-plane ferro— and 
antiferromagnetic^ systems with an in-plane symmetry- 
breaking magnetic field (or anisotropy). Quasi-ID com- 
pounds, obtained by the molecular synthetic approach^, 
now allow the experimental investigation of dynamic 
phenomena that had been theoretically predicted some 
decades ago^. Such systems have been named Single- 
Chain Magnets (SCM's)^'^, by analogy with Single- 
Molecule Magnets (SMM's)^, since they show magnetic 
hysteresis with no evidence of 3D magnetic ordering but 
rather due to slow dynamics of the magnetization of a 
pure ID character. It was shown^iiii that the relax- 
ation towards thermodynamic equilibrium in these sys- 
tems is driven by thermally activated single-spin fiips, 



described in first approximation by the stochastic theory 
developed by Glauber— for the ferromagnetic ID Ising 
model. The effect of introducing non-magnetic impuri- 
ties in the chain was investigate d^^i^^i^*^ , too. The first 
system showing SCM behavior was the quasi-lD molecu- 
lar compound of formula [Co(hfac)2NITPhOMe]oo, where 
hfac is hexafiuoroacetylacetonate, NITPhOMe is the 
nitronyl-nitroxide radical 2-PhOMe-4,4,5,5-tetramethyl- 
4,5-dihydro-lH-imidazolyl-l-oxil-3-oxidei^. The repeat- 
ing unit of such molecular chains is formed by two differ- 
ent magnetic centers, a M'^^ (hfac) 2 moiety containing a 
metal (in the present case Co^^), and a radical moiety that 
bridges the metals. The properties of this class of mag- 
netic polymers can then be tuned, with different effects, 
by rationally changing either the metal or the radical 
moieties. In the rich field that has stemmed from this first 
observation, several other compounds presenting SCM 
behavior have been synthesizedi^iiiii^'i^i2ai2ii2^i2^i21, and 
evidences of quantum effects affecting the magnetization 
dynamics at low temperatures have also been reported^^. 
A one-dimensional hysteretic behavior of the magnetiza- 
tion was identified also in different non-molecular struc- 
tures, such as atomic Co nanowires of finite length, dec- 
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orating the steps of vicinal Pt(997) surfaces^^. 

In general, to observe Glauber dynamics, two require- 
ments are necessary: (i) a strong Ising-like anisotropy 
and (ii) a very low ratio of interchain/intrachain mag- 
netic exchange interactions, in order for slow dynam- 
ics to be observed above the transition temperature to 
3D magnetic ordering. In the framework of molecular 
engineering, quasi-lD chains in which the metal center 
is a rare-earth (RE) ion are very appeahng candidates 
to observe a large variety of phenomena, as RE's are 
characterized by very different magnetic anisotropies^l. 
One limiting case is constituted by molecular-based 
quasi-lD compounds of formula [(Gd)(hfac)3NITR]oo 
(where NITR is 2-R-4,4,5,5-tetramethyl-4,5-dihydro-lH- 
imidazolyl-l-oxil-3-oxide, and R is ethyl (Et), isopropyl 
(iPr), methyl (Me), or phenyl (Ph)), where the Gd(III) 
ions are magnetically isotropic. In these systems heli- 
magnetic behavior could be evidenced, due to competing 
magnetic interactions between nearest neighboring (NN) 
and next-nearest neighboring (NNN) spin sites. When R 
= Et, Villain's conjectured^ of a two-step magnetic or- 
dering to a low temperature helical 3D phase, through 
an intermediate chiral spin liquid 3D phase, was con- 
firmed (the two transitions occurring at 1.88 K and 2.19 
K, respectively)^. When passing from isotropic Gd to 
Dy, which presents a strong Ising-like anisotropy, a tran- 
sition to 3D magnetic order occurring at 4.4 K can be 
observed in [Dy(hfac)3NITEt]o o^°'^^ . However, using a 
bulkier NIT(C6H40Ph) radical^ increases the distance 
between chains, so that SCM behaviour is observed with 
no evidence of any 3D magnetic order down to the lowest 
investigated temperature (1.7 K). The two requirements 
(i) and (ii), necessary to observe Glauber dynamics, are 
thus fulfilled in [Dy(hfac)3NIT(C6H40Ph)]cx,. Consis- 
tently, the relaxation time r was found^ to follow an 

A 

Arrhenius law r = tqc'^b^. In order to fully under- 
stand the rich dynamics of these quasi-lD system o^^i^^ , 
in particular the role played by NNN interactions, a 
modelization of the static properties is necessary. In 
this paper, an accurate experimental study of the an- 
gular dependence of the single-crystal magnetization 
of [Dy(hfac)3NIT(C6H40Ph)]oo, combined with ab ini- 
tio estimation of the magnetic anisotropy of the Dy^+ 
ions, is presented. The static magnetic properties of 
the chain system are simulated by means of a classical 
transfer matrix calculation, formulated for a ID Heisen- 
berg model with antiferromagnetic exchange coupling 
and non-coUinear, local anisotropy axes. Fine agreement 
with experimental data as a function of crystal orien- 
tation, temperature and magnetic field is obtained and 
the Hamiltonian parameters can then be estimated. The 
dominant role of NNN magnetic interactions of these 
molecular structures, despite the large Dy-Dy separation, 
is clearly demonstrated. These findings are expected to 
be extremely important to target novel SCM dynamics in 
molecular compounds, and to rationalize the SCM behav- 
ior of [Dy(hfac)3NIT(C6H40Ph)]oo and the related class 
of compounds^^. In particular, as the role of natural and 



induced defects is substantially different when competing 
NN and NNN interactions are present, this situation can 
lead to the observation of different dynamic regimes in 
addition to the already known ones. 



II. CRYSTAL STRUCTURE AND AB INITIO 
CALCULATIONS 

Recently we investigated in detail a monomeric com- 
pound of formula Dy(hfac)3(NITR)2 where the metal 
ion is coordinated to the oxygen atoms of two different 
radicals^, like in the chain. In fact, we suggest that the 
monomer can be considered as a building block of the 
chain. Using angle-resolved magnetometry, the molec- 
ular magnetic susceptibility tensor of the monomeric 
derivative has been reconstructed, revealing an excep- 
tionally large Ising-type magnetic anisotropy of the Dy^+ 
ion^. This technique also provides the orientation of 
the principal magnetic axes of the monomer; the Ising 
axis turns out to He in between the two radicals, almost 
along the binary axis of the approximate C2 symmetry of 
the molecule. Moreover, the experimental findings nicely 
agree with ab initio calculations within ca. 7° deviation 
between the experimental and theoretical directions of 
the magnetization easy-axis. 

Here, we apply the same quantum chemistry 
method(CASSCF/RASSI-SO^), using the MOLCAS- 
7.0 package^, in order to obtain the single-ion 
easy-anisotropy axes of the Dy^+ ions in the 
[Dy(hfac)3NIT(C6H40Ph)]oo chain compound. In 
this multi-configurational approach relativistic effects 
are treated in two steps, both based on the Douglas- 
Kroll Hamiltonian^. Scalar terms are included in the 
basis set generation and are used to determine spin-free 
wavefunctions and energies, through the use of Complete 
Active Space Self Consistent Field (CASSCF) method. 
Then, Spin-Orbit coupHng is treated with a Restricted 
Active Space State Interaction computation (RASSI- 
SO), which uses the CASSCF wavefunctions as the basis 
states. By using the resulting eigenstates of the previous 
method, the gyromagnetic tensor of the ground doublet 
Kramer's state can be computed and diagonalized in 
order to obtain the three main anisotropy axes and the 
gyromagnetic values along those axes {gx, Qy and Qz). 

Computations are performed in a quantum cluster 
model considering a Dy^+ ion and its surrounding ligand 
molecules, i.e. three hfac and two NIT radicals, using 
the geometry determined by the X-ray structure analy- 
sis. In this cluster model several modifications have been 
introduced in order to reduce the large computational 
time: the fluorides in the hfac ligands have been replaced 
by H atoms, the NIT(C6H40Ph) radicals have been re- 
placed by NIT(C6H5), and H atoms have been added to 
the external oxygen atoms of the radicals transforming 
them in close-shell molecules, and therefore reducing the 
active space of the CASSCF calculations. Despite the 
previous modelization, the quantum cluster would still 
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Table I: Ab mi£jo-calculated energies of the eight doublets 
of the ^Hi5/2 multiplet of the Dy^''' ions in the coordination 
environment determined from the crystallographic structure. 



doublet 


EO 


El 


E2 


E3 


E4 


E5 


E6 


E7 


Energy (K) 





58.3 


84.6 


125.9 


164.9 


238.4 


328.9 


591.0 



require a large computational time due to its size. There- 
fore a Cholesky decomposition technique^, recently im- 
plemented in the MOLCAS package, is used to overcome 
problems arising from the large size of the two-electrons 
integral matrix. 

All the atoms were represented by basis sets of atomic 
natural orbitals from the ANO-RCC library as imple- 
mented in the MOLCAS-7.0 quantum chemistry package 
and the following contractions were used: [8s7p4d3f2g] 
for Dy, [3s2p] for 0, C and N, and [2s] for H. The 
CASSCF active space consisted on the Dy 4f orbitals, 
containing 9 electrons in 7 orbitals. The previous ab ini- 
tio computation in the monomeric compound showed a 
minor role of the quadruplets and the doublets CASSCF 
states in the Spin Orbit coupling effect on the energy 
levels of the ^Hi5/2 ground multiplet. Therefore, only 
the three sextets, ^H, ®F and ®P, have been computed in 
three different CASSCF state average calculations, one 
for each of the sextets. The 21 resulting CASSCF states 
were introduced in a RASSI-SO state interaction in order 
to compute the effect of the Spin-Orbit coupling. The 
computed energies of the ^Hi5/2 multiplet are listed in 
Table H 

The gyromagnetic factors for the ground doublet have 
also been evaluated and found to be gx=0-S, gy =0.7, and 
gz=l&.7. The computed single-ion anisotropy axis al- 
most coincides with a binary axis in the idealized sym- 
metry of the DyOs polyhedron, similarly to what ob- 
served for the monomer. The easy axis forms an an- 
gle 9 of ca. 80° with respect to the b crystallographic 
axis (along the chain axis), whereas the projection of the 
single-ion anisotropy axis in the ac plane (perpendicular 
to the chain axis) forms an angle (j) of ca. 50° with the a 
axis. Compared to the triclinic monomer, here the situa- 
tion is complicated by the fact that the chain crystallizes 
in the P2i2i2i space group (N°19), Z= 4, with three two- 
fold screw axes as symmetry elements. The asymmetric 
unit contains one Dy^+ ion, and consequently all dyspro- 
sium atoms are symmetry related. When superimposing 
the Ising axis of the Dy^+ ion on the chain, its symmetry 
elements generate a canted structure as depicted in Fig. 
[1^. In Fig. [TJd, a view along the chain direction (the crys- 
tallographic b axis) of the crystal packing of the chains is 
reported. As a consequence of this crystal packing, even 
if the chains are structurally all equivalent, they are dif- 
ferently oriented from the magnetic point of view. Two 
types of chains (A and B) are related by a 2i two-fold 
screw axis (in the ac plane), while each of these chains 
is generated by its own 2i axis (along b). When analyz- 
ing the magnetic behavior of the chain, we thus have to 




Figure 1: (color online) a) Representation of the Ising axes 
of the Dy'^^ ion, as gathered from ab initio calculations and 
superimposed on the chain in the ab plane, b) View along 
the b crystallographic axis of the crystal packing, showing the 
two symmetry related types of chains. The orientation of the 
X and Z axes in the ac plane is also represented (the Y axis 
coincides with 6). 

consider not only the canted structure generated by the 
alternate incHnation of the Ising axes with respect to the 
chain direction, b, but also the two families of chains (A 
and B), whose projections on the ac plane of the easy 
axis of the Dy^+ ions form an angle of ca. 80° . 



III. MAGNETIC MEASUREMENTS 

Single-crystal magnetic measurements were performed 
using a homemade horizontal rotator that allows, after 
orientation of the crystal, to measure the magnetic sus- 
ceptibility along and perpendicularly to the chain axis. 
The morphology of the crystals, with large (101) and 
(101) faces, does not allow to perform rotations along 
the principal crystallographic axes. We thus define a 
laboratory frame, {X, Y, Z), with Z corresponding to the 
normal to the (101) crystal face, Y = b, and X orthog- 
onal to the first two. The new reference frame is thus 
obtained performing a rotation of the crystallographic 
frame, (abc), around b (the chain axis), by an angle 
a « 49.8° (see Fig. [lb). 

In Fig. [21 we show the angular dependence of the ra- 
tio M/H between magnetization and field, from here on 
indicated as the molar susceptibility xm, in three orthog- 
onal rotations performed at 2.8 K in an external field of 1 
kOe. As expected for an orthorhombic system, the crys- 
tallographic axes correspond to relative extrema. Two, 
out of three rotations, are almost identical, and display 
a strong angular dependence of xm , which has its max- 
imum along the chain direction. The third rotation, the 
one around 6, shows a very low, angle-independent value 



4 




90 135 180 225 270 315 

Angle in 



Figure 2: (color online) Angular variation of the molar suscep- 
tibility, xm, for three orthogonal rotations at 2.8 K, measured 
in a 1 kOe external field. Rotations were performed around 
X (green circles), b (blue triangles) and Z (red squares). 



of XM- 

The temperature dependent magnetic behavior of 
[Dy(hfac)3NIT(C6H40Ph)]oo was also investigated. 
Measurements were performed in the 2.4-51 K tempera- 
ture range, both for the angular variation of xm (Fig. [3]) 
and for the field (H) dependence of the magnetization M 
(Fig. H]). If we consider the rotation around X reported 
in Fig. mtop) for temperatures between 2.4 and 6.5 K, 
the maximum of xm is found at 90°, thus along the chain, 
and the minima at 0° and 180°, i.e. perpendicularly to 
the chain. At 8 K, the xm curve has almost no angu- 
lar dependence. For temperatures from 9.5 K to 51 K, 
the maxima and minima are inverted with respect to the 
low temperature curves. This is a quite spectacular and 
peculiar feature which, at first glance, seems troubling. 
However, previous investigations^^^ give us some hints 
to rationaHze this result: 

• The direction of the Ising anisotropy of the Dy^+ 
ion in the monomer— is temperature independent 
in the temperature range investigated here, thus 
the anisotropy crossover observed around 8 K can- 
not be attributed solely to an electronic effect of 
the Dy^+ ions in the chain compound. 

• The projection of the Ising axes, gathered from 
ab initio calculations and the analysis of the 
monomeric building block^ of the chain structure 
suggests that the easy axes are not collinear. 

• The observed behavior (i.e., the anisotropy 
inversion) recalls the one encountered in 
[[Mn(TPP)02PPh]-H20]oo (where TPP=meso- 
tetraphenylporphyrin), a Mn^^^-based canted 
antiferromagnetic SCM^. 

It is interesting to note that similar inversions (at two 
different temperatures) were observed^ in the magnetic 
torque curves of TPP[Fe(Pc)(CN)2]2, a molecular con- 
ductor that can be described by a one dimensional 
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Figure 3: (color online) (Top) Angular variation of the molar 
susceptibility, xm ~ ^'^/H, for temperature T ranging from 
2.4 K (blue) to 51 K (red), measured in an external magnetic 
field of 1 kOe. Rotations were performed around the X axis, 
with 0° and 180° corresponding to the field aligned perpendic- 
ularly to the chain along -Z and Z, respectively, whereas, at 
90°, Xm is measured along the chain. Temperature color map- 
ping is described on the right part of the figure. Vertical dot 
lines are guides to the eyes. (Bottom) Transfer matrix simu- 
lation of the angle dependence of XAf = M/H (with H = 1 
kOe) when rotating around X. The classical spin Hamilto- 
nian, Eq. ([T]), and the Hamiltonian parameters specified in 
Section IV were used for the calculation. Lines of different 
colors refer to different temperatures. 



anisotropic Heisenberg model with antiferromagnetic ex- 
change interactions. 

A scenario of a canted AF ID structure similar to the 
one observed in the [Mn(TPP)02PPh]-H20 chain^i re- 
quires that, in the alternating Dy-radical chain under 
study, the NNN AF interactions dominate the NN ones. 
To confirm this hypothesis, the field dependence of the 
magnetization was measured. In Fig. [DJtop) we report 
the M vs H curves recorded along and perpendicularly to 
the chain. A very different shape of the curve is observed 
depending on whether the field is ahgned along the non- 
compensated moment (i.e., along the chain), thus giv- 
ing a rapid saturation of the magnetization, or closer to 
the compensated easy axes (i.e., perpendicularly to the 
chain) , providing a sort of metamagnetic transition when 
the magnetic field overwhelms the antiferromagnetic in- 
teraction (see Section [IV] for quantitative details). 

The rapid saturation of the component along the chain 
observed at 3 K, with a plateau corresponding to ca 10000 
emu mol~^ (« 1.8^^), is the typical behavior expected 
for a weak ferromagnet (WF) - indeed in this case a 
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Figure 4: (color online) (Top) Field dependence of the mag- 
netization along the chain (full squares) and perpendicular to 
the chain (along the Z axis, open circles), for temperature 
T ranging from 2.45 K (blue) to 12.85 K (red). The tem- 
perature color mapping is depicted on the right. (Bottom) 
Transfer matrix simulation of the field dependence of M vs H 
along the chain direction (dashed lines) and perpendicular to 
it (solid lines) obtained using the classical spin Hamiltonian, 
Eq. ([1]), and the Hamiltonian parameters specified in Section 
IV. Lines of different colors refer to different temperatures. 



ID WF - along the direction of the non-compensated 
moments. Perpendicularly to the chain, and more pre- 
cisely along the Z axis, the magnetization is almost linear 
at low field, but shows a rapid increase around 15 kOe 
and saturation at larger fields, almost 22000 emu mol^^ 

(« 4.1^5). 

Interestingly, the step in the magnetization perpendic- 
ular to the chain becomes smoother when temperature is 
increased, and finally disappears for temperatures above 
10 K, when the weak AF interaction becomes neghgi- 
ble and the anisotropy of the sample is driven by the 
single-ion anisotropy of the Dy^+. When the field is 
applied along the chain direction, i.e. the direction of 
the non-compensated moment, a more regular behavior 
is observed, as the initial susceptibility monotonously de- 
creases with increasing temperatures. 

This feature sheds some light on the quite unique 
behavior of the whole family of lanthanide-based 
SCM's^i^^. Magnetization curves recorded on the 4/- 
based SCM family always display steps at relatively high 
field, and a trend was observed along the lanthanide 
series^i^. In fact, the position of the step is explained by 
the simultaneous presence of an anisotropic center (pro- 
viding the anisotropic canted structure) and of an AF 



interaction along the chain. Hence, the stronger is the 
intrachain magnetic interaction, the higher is the field 
needed to compensate this interaction, shifting the step 
towards higher fields. 

This allows concluding that, in the 
[Dy(hfac)3NIT(C6H40Ph)]oo SCM under study: 

• Antiferromagnetic interaction between Dy^+ ions is 
present. Due to this fact, the magnetization com- 
ponents perpendicular to the chain direction cancel 
out at low temperature and low field, while they 
become favored at higher temperatures. 

• The high temperature hard axis becomes the easy 
axis at low temperature because non-compensation 
of the canted spins occurs. 

• The anisotropy inversion is visible around 9 K, sug- 
gesting a relatively weak intrachain exchange inter- 
actions, as expected in compounds involving tripos- 
itive lanthanides. 

If we now compare the observed magnetization values for 
the chain and the monomer derivative, a first saturation 
-^'-^sat ~ 58000 emu mol^^ (« II^b) has been observed 
along the easy-axis of the monomer—. At a first level of 
approximation we can neglect the radical contribution, 
which is much smaller than that given by the paramag- 
netic ions; moreover, the "radical" sublattice is on itself 
experiencing AF interactions mediated by the lanthanide 
ion and is also frustrated due to the sizeable AF NNN in- 
teraction between Dy^+ spins (see Section [TV|> ^. This al- 
lows providing a rough estimation of the angle, 9, formed 
by the easy axis of each Dy^+ center with the chain axis, 
b. In fact, comparing the afore-mentioned saturation 
value (« 58000 emu mol~^) with the rapid saturation one 
(w 8000 emu-mol~^) reached when the field is applied 
along the chain, yields 9 = cos-^ (8000/58000) « 82°. 
This value is consistent with the results of the ah initio 
calculations. 

As for the magnetization curve measured perpendicu- 
larly to the chain, the field was appHed aligned along Z. 
Considering the direction of the local easy axes predicted 
by the ah initio calculations (see Fig. [TJd), the applied 
field is almost parallel to the projection of these axes 
on the ac plane for type B chains, whereas it is almost 
perpendicular to type A chains. For the two symmetry- 
related families of chains, A and B, the projections of the 
local easy axes on the ac plane are almost reciprocally 
perpendicular but, remarkably, this is a fully acciden- 
tal coincidence rather than a symmetry-imposed result. 
Therefore, as schematically represented in Fig. [H the 
applied field can compete with the AF interaction inside 
the chain of type B whereas it has no effect on the other 
family (A), where Dy^+ spins remain antiferromagneti- 
cally coupled. Thus, only one chain type is contribut- 
ing to the magnetization, which should approach a value 
M = (i)58000sin6l (that is, M « 28700 emu-mor^), not 
very far from the experimental value. Despite the fact 
that these figures can be only taken as indicative, given 
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Figure 5: (color online) Influence of an applied magnetic fleld 
in the ac plane on the orientation of the magnetic moments 
of neighboring Dy ions of the two types of chains (A and B) . 



the rough approximation of neglecting the radical contri- 
bution, such saturation values for the magnetization are 
in agreement with our prediction about the orientation 
of the local easy axes, based on the ab initio calculations. 

A further confirmation of the correctness of this hy- 
pothesis can be afforded by performing, at low tempera- 
ture, a rotation around b (the chain axis) under an exter- 
nal field of 30 kOe: i.e., a field strong enough to be capa- 
ble of distinguishing the different projections of the Ising 
axes in the ac plane by overcoming the AF interaction. 
From Fig. (Utop) one can clearly see that, contrarily to 
what observed at low field, xm is now angle dependent, 
with a very peculiar behavior with large maxima and a 
roughly 90° periodicity. It should be noticed that: 



• The absolute maxima are observed along the crys- 
tallographic axes a and c, and not along the direc- 
tions where the projections of the local axes fall. 



• The data show a secondary periodical structure of 
smaller maxima shifted by an angle of ca. 45°. 

As a final remark, we notice that the maxima observed in 
Fig. 6 (top) along a and c are not strictly equivalent, but 
the direction c appears to be favored. In the following 
Section we will show that this feature can be interpreted 
as due to a slight (but appreciable) departure from or- 
thogonality of the projection on the ac plane of the local 
easy axis of Dy^"*" ions belonging to the two symmetry- 
related chain families, A and B, as indeed predicted by 
ab initio calculations. 



We have also checked the temperature dependence of 
the anisotropy when rotating at high field around the 
chain axis, b (see Fig. [71 top). The results show that the 
behavior is strongly temperature dependent. Since in this 
temperature range the Dy^"^ anisotropy of the monomeric 
constituent units of the chain does not change, these data 
confirm the key role played by the exchange interaction, 
suggesting an order of magnitude of the AF interaction 
of a few Kelvins (see the following Section for a more 
detailed discussion). 
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Figure 6: (color online) (Top) Angular variation of the molar 
susceptibility (xm ~ M/ H) for the rotations around the chain 
axis h at 2.5 K, measured in a 30 kOe external fleld (blue tri- 
angles) or in a 1 kOe external fleld (black triangles). (Center) 
Transfer matrix simulation of M/H vs Angle at 2.5 K in a 
fleld of 30 kOe, showing the contribution from each chain type 
(A and B), as well as their sum. The classical spin Hamilto- 
nian, Eq. ([T|), and the Hamiltonian parameters specifled in 
Section IV were used for the calculation. Perfect orthogonal- 
ity 2(/)=90° between the projections on the ac plane of the 
easy axes of the two chain families A and B was assumed. 
(Bottom) Schematic representation of the corresponding spin 
conflgurations. 



IV. THEORY AND DISCUSSION 

In order to interpret the measurements of the static 
magnetic properties of the chain, we model both A and 
B chains with a classical spin Hamiltonian 

n = -^{JS(z)-S(* + l) + i?[5,.(z)]2 



a,i3 



(1) 



7 




90 135 
Angle / (°) 



Figure 7: (color online) (Top) Angular variation of the sus- 
ceptibility XM = M/H from 2.45 K (blue) to 20.45 K (red) 
measured in a 30 kOe external field. Rotations were per- 
formed around the chain axis b as defined in the text. Color 
mapping is described on the right part of the figure. Solid lines 
are guides to the eyes. (Bottom) Transfer matrix simulation 
of XM = M/H vs Angle, at selected temperatures in a field 
of 30 kOe. The classical spin Hamiltonian, Eq. (HJ, and the 
Hamiltonian parameters specified in Section IV were used for 
the calculation. Notice that, with respect to Fig. [Sljcenter), 
a slight deviation from perfect orthogonality between the two 
chain families was assumed, thus accounting for asymmetry 
in the height of the major peaks. 



where is the Bohr magneton and a, /3 = X,Y, Z de- 
note Cartesian components in the crystallographic frame 
{a,b and c), while the subscript z- denotes the spin com- 
ponent along a local anisotropy axis, varying with lattice 
site i. In the sum, only the Dy centers are considered 
for the reasons given below. The AF exchange interac- 
tion (J < 0) is considered isotropic, while the anisotropy 
is introduced through a single-ion term {D > 0) favor- 
ing local z[ axes, which are different for even and odd 
sites, and in the Lande tensor (ga/s)- In a local reference 
frame, the Lande tensor is assumed to have only diagonal 
components, and g± (where parallel and perpendicu- 
lar refer to the local axis z'J. Since the spin Hamilto- 
nian U) is classical, the spin modulus S does not matter 
(its effect can be taken into account by the parameters 
J, D, and g's). For an easy comparison with the cal- 
culation performed for the monomei>^, we will assume 
S = 1/2. Notice that the classical nature of the spins 
allows the inclusion of single-ion anisotropy even in this 
case^. Concerning the neglect of the radicals, we note 



that their anisotropy is very low compared to the Dy^+ 
one, and the Dy-radical interaction is necessarily frus- 
trated, since the interaction between the Dy^+ ions is 
antiferromagnetic. As a first approximation, the system 
can thus be considered as the superposition of two inde- 
pendent one-dimensional magnetic lattices, one formed 
by the Dy^+ ions (with magnetic anisotropy) and the 
other by the radicals (magnetically isotropic). At low 
temperatures, the contribution of the "radical" lattice is 
expected to be very small because of the relatively strong 
NNN antiferromagnetic interaction between radicals^. 

Using a transfer matrix method^, the static magnetic 
properties of the ID classical spin model |(T|) can be ex- 
actly expressed in terms of the eigenvalues and eigen- 
vectors of a real, positive-definite matrix (see Appendix 
A for details). In this way, starting from certain values 
of the Hamiltonian parameters, we were able to numer- 
ically simulate the static magnetic behavior of the real 
chain compound. The simulation procedure started from 
placing the local easy axes at 6* = 75° from the b axis and 
its projection on the ac plane aX (j> = ±45° from a. The 
other Hamiltonian parameters used for the simulation 
were J = -24 K, D = 160 K, .g,, = 21 and g± = 4, where 
parallel and perpendicular refer to the local z- axis. The 
contributions from the two families of chains were then 
summed up. In Fig. [3l[bottom) we report a simulation of 
the quantity M/H, calculated for different temperatures 
when rotating along AT, to be compared with the mea- 
sured dc susceptibility in Fig. [3l[top). The susceptibility 
presents a maximum along b and a minimum along Z for 
T < 8 K, while the trend is reversed for T > 8 K. Thus, 
the anisotropy inversion - observed at T = 8.5 K on the 
experimental data in Fig. [Sjtop) - is nicely reproduced. 

The field dependence of the magnetization at selected 
temperatures, both for H parallel and perpendicular to 
the chain, was simulated using the same Hamiltonian pa- 
rameters: see Fig. IDJbottom). The measured overall be- 
havior of Fig. mtop) was reproduced, although a more 
quantitative comparison was prevented by the neglect of 
the radical contribution in the adopted model U). 

In particular, a relationship between the field at 
which a rapid increase of the magnetization perpendic- 
ular to the chain is observed and the Transfer-Matrix- 
computation parameters can be deduced by simple argu- 
ments. Without loss of generality, we assume that spins 
are essentially constrained to be oriented along their lo- 
cal anisotropy axes at low temperatures and consider an 
applied field lying on the ac plane at = -1-45° from a. 
This implies that the Zeeman energy of half of the spins 
is frustrated (the S^'^^-spin sublattice with the notation of 
Appendix A). The fiip of such spins lowers the Zeeman 
energy per spin by a factor AEzeeman = 2Sg^^^BH svaO . 
On the other hand, the corresponding increase of the 
exchange energy per spin is AEexch = 4S'^| Jcos(20)|. 
The reversal of such a sublattice is expected to occur 
when these two energies are equal, i.e. at the field 



H 



flip 



2S|,/cos(2fl) 
g.. fiB sin 6 



15.25 kOe. These naive energy 
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considerations give the correct value for the "jump" field 
(see Fig. |4lJbottom)), the effect of temperature being 
solely that of rounding the step in the magnetization 
curve. 

Under quite general conditions^, the correlation 
length £, of anisotropic, classical-spin chains behaves 

asymptotically at low temperatures as ^ ^ exp (^-^f^j^^ • 

Consequently, for T > both famihes of spin 

chains are expected to lose their genuine ID charac- 
ter and respond to a field as independent, anisotropic 
paramagnetic centers. Again, the order of magnitude of 
^2kl'''' ~ ^^-^ K matches with the temperature at which 
the anisotropy inversion was observed. 

The rotation around the chain axis, b, of the molar sus- 
ceptibility was also simulated in a strong external field 
{H = 30 kOe), and the contributions from each fam- 
ily of chains (A and B) are reported in Fig. [6lJ center), 
along with their sum. A qualitative explanation of the 
observed behavior can be obtained by looking in detail 
at the effect of the field when it is applied at different 
angles in the ac plane: see Fig. [Bljbottom). When the 
field is along a (or c), it lies at nearly 45° from the easy 
axis of both famihes of chains. Even if a smaller compo- 
nent of the field (indeed proportional to cos 45°) acts on 
the magnetic moments, this is still enough to overcome 
the AF interaction. This is for example represented by 
the arrows at Angle = 45° in Fig. [^bottom). Given the 
large anisotropy of Dy^""" in this environment, the mag- 
netic moments remain aligned along their easy axis, and 
only the weak intrachain AF interaction is overcome. The 
largest magnetization value we can measure is therefore 
reduced to Msat sin 75° cos 45° « 39600 emu mol~^, not 
too far from the observed value {~ 37500 emu mol~^). 
It is evident that each chain type gives a contribution 
with a periodicity of 180° so that, owing to the acciden- 
tal orthogonality between the two families A and B, the 
overall periodicity is « 90°. 

Finally, in order to simulate the asymmetry in the 
height of the major peaks displayed in Fig. [SJtop), the 
hypothesis of a perfect orthogonality between the projec- 
tions of the two chain types, A and B, was released. In 
Fig. [T^bottom) the angular variation of xm at selected 
temperatures, calculated assuming an angle (j) between a 
and the projection of the local easy axes in the ac plane 
for the A-type chain slightly larger than 45°, is reported. 
The fine agreement with experimental data (see Fig. [3 
top) proves that the different height of the major peaks 
at 45° and 135° can be attributed to non-perfect orthog- 
onality. Interestingly, even the minor maxima between 
two main maxima are reproduced at sufficiently low tem- 
peratures. 



V. CONCLUSIONS 

In conclusion, we have performed a thorough ex- 
perimental and theoretical investigation to rational- 



ize the complex magnetic behavior of the Dy-based 
SCM [Dy(hfac)3NIT(C6H40Ph)]oo. By means of angle- 
resolved magnetometry the magnetic anisotropy of this 
complex system has been accurately determined and it 
has been possible to evidence a peculiar "anisotropy- 
inversion" phenomenon, in which the easy and hard axes 
of the magnetization of the system swap when rising the 
temperature. These features could be fully explained, 
with excellent agreement between theory and experi- 
ments, using a combined theoretical approach in which 
the single RE-ion anisotropy is described at the quantum- 
chemical level while the thermodynamic properties of the 
whole system of coupled ions are computed with a classi- 
cal spin Hamiltonian. To our knowledge, this is the first 
report of such a combined theoretical strategy, which we 
found indeed mandatory to tackle the complexity of some 
magnetic molecular materials, often characterized by a 
low symmetry environment of the magnetic centers. The 
rich magnetic behavior of [Dy(hfac)3NIT(C6H40Ph)]oo 
was well reproduced by the model calculations that high- 
lighted the presence of two sets of non-interacting par- 
allel chains, with their projections mutually tilted by 
about 90°. Accordingly, the metamagnetic transition 
inside each symmetry-related family of chains could be 
addressed independently. The combined approach has 
also allowed quantifying the sign, magnitude and role of 
NNN interactions between Dy'^+ ions. Despite the large 
distances between the RE's (> 8 A), the overall behav- 
ior is actually dominated by this antiferromagnetic NNN 
exchange pathway, which gives rise to the very rich mag- 
netic behavior observed. The same evidence was previ- 
ously reported for isotropic Gd ions. The present find- 
ing shows that NNN are more common than believed so 
far and may induce a critical revision of the description 
of the magnetic properties of polynuclear compounds. 
The present study may thus represent a methodological 
guide fine to rationalize the magnetic properties of RE- 
based SCM's as well as other extremely complex behav- 
iors displayed by molecular materials in general. To the 
aim of understanding the dynamic behavior of RE-based 
SCM's, the obtained information represents a fundamen- 
tal, preliminary step; from the knowledge that compet- 
ing NN and NNN interactions are present, for instance, 
one can reasonably expect the observation of novel dy- 
namic regimes. In particular, we expect the effect of 
natural and induced defects to be substantially different 
than the already known cases in which only NN inter- 
actions are presen t^^i^^i^"^ . As an example, preliminary 
results indicate that, while the contribution of the rad- 
icals to the static properties is almost negligible, they 
actively contribute to the correlation length even when 
the Dy'^+ ions are partially substituted by diamagnetic 
ions. In this sense, by doping with diamagnetic ions the 
controlled interaction between segments of chains can be 
tuned, thus allowing a further form of engineering of the 
dynamic properties. 
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I A-type chain 

Figure 8: (color online) Canted ground state of the spins of 
a single chain of type A, described by Hamiltonian For 
each unit cell I, there are two spins (thick, green arrows, rep- 
resenting Dy^+ ions) that belong to two different sublattices, 
/ and //, with different local axes Zj and z'jj (thin, dashed 
arrows). 
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The Hamiltonian of the chain can be rewritten, as in 
Ref. [U, in the form 

N 

H ^ $^{V^i[S^(0, S"(0] + V2[S''{1), S'{1 + 1)]} (A.l) 

1=1 

where the total number of spins, in the single A-type 
chain considered, is 2N because there are two spins per 
each unit cell {I = 1, • • • , N) . In our case, Vi and V2 
take the forms (a, (3 = X,Y, Z) 

Vi[S'{l),S'\l)]^-JS\l)-S'\l) 
-^Mo E HAgipSj^il) + gi'f,S'/il)} (A.2) 

a, (3 
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Appendix: TRANSFER MATRIX METHOD 

In this Appendix, we briefly describe the transfer ma- 
trix method, which allows the determination, at any 
nonzero temperature, of the static properties of a one- 
dimensional model of classical spins in terms of the eigen- 
values and eigenfunctions of a real, positive-deflnite ma- 
trix. Here the method will be appHed in the thermody- 
namic limit, but it can be reformulated also to treat the 
case of a finite system with open or periodic boundary 
conditions. 

It is worth noticing that the method is exact: the only 
approximation is that the spins are considered as classical 
vectors. Since the method requires an integration over 
the unit sphere surface to be discretized, numerical errors 
are significant only at very low temperatures, where the 
integrand presents strong variations over the integration 
domain. However, the accuracy of the method can be 
improved at will by using a greater number of integration 
points. 

Following Pandit and Tannous^, we consider a single 
magnetic chain with two sublattices, I and // (not to 
be confused with the two families of chains related by 
"accidental" orthogonality, A and B). 

To fix ideas, let us consider a single A-type chain: in 
the l-th unit cell, there are two spins, S^(^) and S^^{1), 
characterized by two different local axes, z'j and z'jj, 
forming equal angles (6*/ = On = 75°, see Fig. [8]) with 
the chain axis b and their projections in the ac plane dif- 
fering by 180° (0/ 45°, = 225°) with the a axis. 



V2[S''{1), S'{1 + 1)] = -JS"{1) ■ S\l + 1) 
-^^0 E H^ig'JpS'/il) + gipSlil + 1)} (A.3) 

a, 13 



Clearly, in the previous equations one has to take 
into account that the crystallographic (laboratory) frame 
(AT, y, Z) and the local frame (x',y',z') are related by a 
rotation, so that for the /-sublattice one has 



Si' 


~ sin0/cos(/)/ S\q + sin 0/ sin (/)/ Sy 






+ cose iSz 


(A.4) 


9xx 


= gx(cos^ 0/ cos^ 0/ -1- sin^ 0/) 






+ sin^ 9 1 cos^ (pi 


(A.5) 


9XY 


~ (ffll ~ sin^ 9 1 cos (jyj sin 0/ 


(A.6) 


9xz 


~ (.9|| ^ cos 9 1 sin 6'/ cos </>/ 


(A.7) 


9yx 


= 9XY 


(A.8) 


9yy 


~ g± (cos^ 9 1 sin^ (j)i + cos^ (/>/) 






+ sin^ 9 J sin^ 


(A.9) 


9yz 


~ {g\\ — g_L) cos 9i sm9i siiKpi 


(A.IO) 


9zx 


9xz 


(A.ll) 


9ZY 


= 9xz 


(A.12) 


9zz 


= g± sin^ 9 1 + g\\ cos^ 9i 


(A.13) 



and similar equations for the other (//) sublattice. 
The partition function Z = Tr e^^"^ can be written as 

Z = j dS^(l) j dS^\\)--- j dS\N) J dS^\N) 



X e 
X e 



-/3yi[S^(l),S"(l)]g-;9V2[S"(l),S^(2)] . . , 
-/3Vi[S^(Af),S"(JV)]g-/3V2[S"(iV),S^(l)] 



(A.14) 
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where periodic boundary conditions, S^(A'' + 1) = S^(l), 
were assumed^. A symmetric, positive-definite l-cernel T 
can now be defined 



(A.15) 



Denoting by A„ and ^'„[S^(Z)] the eigenvalues and corre- 
sponding eigenfunctions of (|A.15|) 

J ds'{i + 1) r[s^(/), s^{i + 1)] *„[s^(z + 1)] 

A„^'„[S^(0] (A.16) 

one has that the largest eigenvalue Ai is always non- 
degenerate, and the eigenfunctions satisfy the two condi- 
tions 

T[S^(0,S^(/ + 1)]=^A„ 

n 

xK[S'mn[SHl + l)] (A.17) 



(A.18) 



of completeness and orthonormality, respectively. Now, 
using first l|A.17P to expand the kernel (jA.lSp . and then 
exploiting (jA.lsp , the partition function Z can be rewrit- 
ten as 

Z=J dS'{l) J dS\2)--- j dS^{N)x 

T[S^(1), S^(2)]r[S^(2), S^(3)] • • • T[S'iN), S'{1)] 

oo 

= ^(A„)^^A[ 



(A.19) 



n=l 



where the largest eigenvalue Ai dominates in the thermo- 
dynamic limit N oo. Using a similar procedure, the 
crystallographic components (a = X, Y, Z) of the mag- 
netization per site on the /-sublattice can be expressed, 
for ^ oo, as 

(Si) = J dS'^US')Si^i{S') (A.20) 

In order to calculate the same quantities on the other 
(//) sublattice of the single (A-type) chain considered. 



one defines another symmetric, positive-definite kernel U 
and another integral equation, respectively 



U[S''{l),S"{l + l)] = J dS'{l + l) 



xe 



-l3V2[S"(l).S'(l+l)] -l3Vi[S'(l+l),S"{l + l)] 



(A.21) 



=A„$„[S^^(/ + l)] (A.22) 

Notice that the eigenvalues of (|A.22p and l|A.16P are 
equal, so that the partition function Z is the same, while 
the eigenfunctions are different. Similarly to Eq. (jA.20p . 
the crystallographic components {a = X, Y, Z) of the 
magnetization per site on the //-sublattice can be ex- 
pressed, in the thermodynamic Hmit, as 

{S'J) = j dS" 'PUS" )S'J^ lis") (A.23) 

For small applied fields, the crystallographic components 
(a = X, Y, Z) of the molar susceptibility xlj for the 
single (A-type) chain considered are then obtained from 
(|A3o1) and (|AT23l) as 



Xm 



^^oA^o;^ E {aUsD + qM')) (A.24) 



where A'o is Avogadro's number. 

The total susceptibility, to be compared with experimen- 
tal data, is then obtained by repeating the calculation for 
a chain belonging to the other family^, i.e. B (related 
to A by "accidental" orthogonality), and then averaging 
over the two contributions. 

As regards the details of numerical calculations, fol- 
lowing Pandit and Tannous^, the integrals over the 
surface of a sphere in (|A.16P and l|A.22|) were approx- 
imated by using McLaren's 72-points, 14th-degree for- 
mula. Thus, the eigenvalues and eigenfunctions of the 
two integral equations l|A.16P and (jA.22p were obtained 
solving a (72 x 72)-matrix eigenvalue problem by standard 
numerical methods^. Because of discretization, errors 
are significant only at very low temperatures, where the 
integrands in (|A.16P and l|A.22p present strong variations 
over the integration domains. However, the accuracy of 
the method can be improved by using a greater number 
of integration points. 
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